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Abstract 

We give a simple and self contained introduction to quaternions and their prac- 
tical usage in dynamics. The rigid body dynamics are presented in full details. 
In the appendix, some more exotic relations are given that allow to write more 
complex models, for instance, the one of a satellite with inertial wheels and 
expressed in a non-inertial reference frame. As it is well known, one nice ad- 
vantage of quaternions over Euler angles, beside the usual arguments, is that it 
allows to write down quite complex dynamics completely by hand. 
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1 Quaternions 
1.1 Fundamentals 

Relation (1), together with associativity and distributivity is all what we will 
use to derive the basic practical applications for quaternions. 



i 2 = j 2 = k 2 = ijk = -1 



(1) 



By left- and right-multiplication in the above equation, we can write 



i ijk = —jk = —i 
ijk k = —ij = —k 
j jk = -k = ji ij j = -i 
i ij = —j = ik ji i = —j 



kj 
—ki 



This shows the product is non commutative and gives the basic multiplication 
rules: 



ij = k 


ji = 


-k 


jk = i 


kj = 


—i 


ki = j 


ik = 


-j 



(2) 



1.2 Notations and Definitions 

A quaternion q is a set of four parameters, a real value qo and three imaginary 
values qii, q 2 j, q^k with qi,q 2l qz E R; it may be written 

q = q + qii + q 2 j + qs,k- 

However, this notation proves itself to be very unpractical. We will therefore 
use two different notations: 

• The quaternion q as a pair of real value and vectorial imaginary value 
q=(q ,q) Rc{q} = q a Im{q} = q = (q 1 q 2 q 3 ) T 

• A column vector of four parameters 

q = (50 qi <?2 <Z3) T 

The conjugate q of q is defined as 

q = (qo,-q) 

and it's norm (a nonnegative real value) as 

\q\ = |q| = \/qi + q 2 i+qi + ql 

The product of two quaternions written as pairs, as described in the next 
section will be noted with o. 
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1.3 Quaternion Product 

From the rules given in (2), we may write the product of q with p. 



(<7o + qii + qi3 + qzk){p + p x i + p 2 j + p 3 k) 



Poqo + 


qoPl i 


+ ?0P2 3 


+ qoP3 k 


+ qipo i + 


q\P\ a 


+ q\P2 ij 


+ q\P3 ik 


+ ?2P0 3 + 


qiP\ ji 


+ 92P2 33 


+ ?2P3 jk 


+ qzPo k + 




+ 93P2 kj 


+ <?3P3 kk 


Poqo 


- qm 


— ?2P2 — 


93P3 


+ {qiPo 


+ qoPi 


+ 92P3 - 


93P2) i 


+ (q2Po 


+ qoP2 


+ <?3Pl ~ 


?iPs) 3 


+ (q3Po 


+ qoP3 


+ qiP2 ~ 


qiPx) k 


qop = 


{poqo -p- 


q, q p + p q + 


qxp). 



(3) 

From (3) it turns out that 

qoq = qoq=(\q\ 2 ,d) = \q\ 2 (4) 
and if q is normed (\q\ = 1) 

qoq = qoq = (1,0) = Id. (5) 

In (3) we also see that 

q o p — p o q (6) 

that is 

\qop\ 2 = (gop)o (fop) = gopopo^ \p\ 2 {qoq) = \q\ 2 \p\ 2 

\P\ 2 

\qop\ = \q\\p\. (7) 

1.4 Quaternions and Spatial Rotations 

First, note the following relations 

(u x v) X w = (u ■ w)v — (v ■ w)u 

9 09 1 — COS <p n ip 1 + cos w 

sin 77 = 7: — - cos 2 £ = — -. 

2 2 2 2 

From now on, q will generally represent a normed quaternion (\q\ = 1) 
involved in a rotation. Let's now place a vector f £ R 3 in the imaginary part 
of a quaternion x and see what happens with it in the following relation 
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x' = qoxoq x = (0,3?) q=(q Q ,q). 

Using (3) 



x = (q- x, qax — q x x) o q 
((q- x)q - (q x- qXx) -q, (q ■ x)q + q (q x - q X x) + (q x - qx x) x q) 

v v / S ^ ' 

Rc{x>} Im{x'} 

Re{x'} = (q- x)q - q (x ■ q) - (q x x) ■ q = 

x = (0, ir'), 
lm{x'} = x" 

= (q- x)q + qlx - q (q X x) + qo(x X q) - (q X x) X q 
— (q- x)q + q$x + 2q (x x q) — (qx x) X q 
= (q- x)q + qlx + 2q (x X q) - (q- q)S + (x ■ q)q 
= 2(q- x)q + q^x + 2q (x x q) - (q- q)x. 



A valid normed quaternion (|g| = yj (q% + q\ + q\ + q\ ) = 1) would be 

q = (qo,q) = (cos ~, sin -n) \n\ = 1. 

In this case, x' becomes 

x =2sm — n • x )n + cos — x + 2 cos — sm — a; x n) — sin —x 
2 K ' 2 2 2 y ' 2 

= (1 — cos tp)(ri ■ x)n + cos <p x + sin tp (xy.fi). 

This last relation is the formula for a rotation by an angle tp around a normed 
axis vector ft, as can be shown with the following figure as follows: 



v 2 = cos f vi + sin ip v 3 
vi = x — (x ■ n)ri 

V3 = Vl X ft 

= (x — (x ■ n)n) x ft 
= (x x n) — (x ■ n) 

o 

=> v 2 = cos f (x — (x ■ n)n) + sin ip (x x ft) 

x = {x ■ n)ft + f>2 

— (x ■ n)ft + cos ip {x — (x ■ ft)ft) + sin ip (x x ft) 
= ( 1 — cos ip) (ft ■ x)ft + cos ip x + sin if (x x ft). 

Moreover 



x = qo x o q 
qox'oq = qoqoxoqoq. 

(1,0) (1,0) 

Thus we have the relations for the rotation and its inverse 



x 



qo x o q 



x 



qo x oq . 



(8) 



1.5 Quaternions and Rotation Velocity 



We will now derive the relation between the rotational velocity vector and the 
quaternion time derivative, x' is any constant vector within the body (rotating) 
reference frame and x is the same vector in the fixed reference frame. As seen 
before, both vectors can be put in relation with 



x = q o x o q x 1 = q o x o q. 
Applying the time derivative to x = (0, x), with x' = (0, x') and x! = 0, we 



get 



x = qox'oq + qox'oq 
x=qoqoxoqoq+qoqoxoqoq 
Id Id 

x = qoqox + xoqoq 



(9) 



and from (3) 



G 



5 ° q = ( go go + q ■ q , -qoq + qoq- 5 x q) 

© 

© = gog'o + 5i 9i + <72<?2 + 5393 = q • q = 
because |q| = 1. That is 

5°5= (0, v) and similarly 9 o g = (0, — v). 

1.5.1 Rotation Velocity in Fixed Reference Frame u> 

From (9) and (10) and using (3) we have 



(10) 



x = (0, u) o x — x o (0, —v) 
x = V y. x — x x v = 2v y. x 

and from (7) 

\x\ = \2v\\x\ => v±x 
If x undergoes a pure rotation, we know that 



x = lu x x and 



thus 



cu = (0,d5) = 2(0, v) = 2qoq 



And right-multiplication by q 



(11) 



uoq~2qoqoq => oj o q = 2q 
Id 



q=\ujoq 



1.5.2 Rotation Velocity in Body Reference Frame u>' 

to' = qo to o q with to = 2q o q 
=>■ Lu' = 2qoqoqoq 
id 



uj' = 2qo q 



And left-multiplication by g 



g o to' — 2 q o q oq = 2q 
id 

q= kou' . 



(12) 



(13) 



(14) 
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1.5.3 Matrix-Product Notation for uj 

From 



lo = 2q o q 



and using (3) 



oq] 


= 2(- 


-4oq-\ 


-9o<f- 


Qi 


Qo 






<12 




9o 


-91 


Q3 


-Q.2 


31 


<7n 



/ 90 \ 

?1 

V 93 y 



w = 2Eq. 

Changing the sign and inverting the cross product allows to make an other 
identification 

lj = -2(-q q + q Q q- qx q) 
uj = -2Eq. 

So the rotation velocity vector in the fixed reference frame can be written as 



uj = 2Eq = -2Eq 



(15) 



And from 



q = -uj o q 



one can similarly find 



uj = (0, uj) => ujq = 



1 



q - 2 I (qoZ + Zxq) J " 2~ E " 



q = \E T uj 



1.5.4 Matrix-Product Notation for uj' 

From 



uj' = 2qo q 



and using (3) 



(16) 
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/ = Im{2g oq] = 2(q q-q q-q x q) 



91 


9o 


93 


-92 


92 


-93 


go 


9i 


93 


92 


~9i 


9o 



9i 

92 

V 93 y 



w' = 2Gq. 

Changing the sign and inverting the cross product allows to make an other 



identification 



uj' = -2(q q- q q- 9x9) 



u>' = -2Gq. 

So the rotation velocity vector in the body reference frame can be written as 



w' = 2Gq = -2Gq 



And from 



(17) 



q = -qoLU 
one can similarly find 



J = (0, uJ') =>• = 



(-9 ' w') ^ 
2 V (ffo^' + qx ui') 



q = \G T & 



-G T u7' 



(18) 



1.5.5 Rotation Matrix i? 

We already have 



uj = 2Ec[ = -2Eq 

q = -# T w 
2 

So we can write 



w' = 2Gq = -2Gq 



w = 2£q 

= 2S(i^ T w) 
v 2 y 



uj' = 2Gq 

= 2G(iG T J') 
= GG T uJ' 



EE T = Id 



GG T = Id 
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And by mixing both sides 



uj' = 2Gq = 2G(-E t uj) = GE T u 



Q = 2Eq = 2E{-G T u') = EG uj' . 

We shall now remember that uj is a vector in the fixed reference frame and 
that uj' is the same vector in the body reference frame, that is Q = Bio' . By 
comparing with the previous two results, we find 



B = EG T 



and 



B- 1 = B T = GE 1 



(19) 



1.5.6 Ep and Gp 

From the identifications made in sections 1.5.3 and 1.5.4, we can see that the 
general meaning the product of E and G with any quaternion p is 



Ep = Im{p o g} 



Gp = Im{<7 op}. 



And from 



it follows 



qoq = qoq=(\q\,0) = (1,0) 



(20) 



Eq = 



Gq = 



1.5.7 One Last Relation 

For any v and due to associativity 



(SW) 

2qoq 



o(0, V) = {—UJ 1 ■ V, U)' X v) 

= 2q o q o v 

= 2{qoqo+q- q, qoq-qoq-q x q) ° v = 2q o (q v -q- v, q v+v Q q+q x v) 



( 



= 2 



qo 



qi 



q-i 



q^_\ 



\ 



q 

G 



qi 

q-2 

T 



-93 

qi 



93 

30 

-9i 



q G 1 
2GG T v = fl'v =J'xt- 



-92 

91 

9o 



v 



I 9o 
9i 

92 
V 93 



-9i 
9o 

93 

-92 

— LO 1 ■ V 

uj' X V 



-92 
-93 

9o 

9i 



92 

-9i 
9o 





< \ 




Vl 








\ vz / 
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Comparing with (17), we conclude that 



Of = 2GG T = -2GG T 



and 



Q,'v 



LU' X V 



(21) 



1.5.8 Relations Summary 



The following table summaries the developed relations, q is always a normed 
quaternion, that is q$ + q\ + q\ + of = 1. 



Quaternion notation 


Matrix notation 


Fixed ref 


Body ref 


Fixed ref 


Body ref 


x = q o x' o q 


x' = q o x o q 


x = Rx' 
R = EG T 


x' = R T x 

R T - R- 1 = GE T 


u> = (0, 65) = 2qoq 


J = (0,3') = 2qoq 


lu = 2Eq = -2Eq 


to' = 2Gq = -2Gq 


q = jLO o q 


q = \q o oj' 


q = \E T u 


q = \G T ui' 






EE T = Id 


GG T = Id 


qoq = qoq = (|g|,0) 


Eq = 6 


Gq = 




(0, u) o (0, v) = 

(-U)' ■ V, Ul' x v) 




0' = 2GG T 

= -2GG T 
Vl'v = to' X V 



(-qi qa -qz 92 \ / -qi <7o <Z3 ~32 \ 

-52 93 9o -9i G = -92 -93 9o 9i 

-93 - 92 91 90 / V -93 92 -9i 9o / 

1.6 Rigid Body Rotational Dynamics 

We now will have a look at the dynamics of a freely rotating rigid body to which 
a momentum T' is applied. Translation of the body will not be discussed (it 
can be decoupled from the dynamics of rotation and is fairly easy). We will 
also consider a potential free system, so that the Lagrangian resumes to the 
rotational kinetic energy only 

L = E rot = X -uj' T Juj' . (22) 

Using the quaternion q as coordinates and with the constraint C — q T q = 1, 
Lagrangian dynamics gives 
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dt dq dq " q ^ dq ' 

F q is the 4-vector of generalized forces which will be expressed in term of 
applied torque later. A is the Lagrangian multiplier used to satisfy the constraint 
C. 

1.6.1 Derivatives of L 

Note the following reminder 

^ = A 
ax 

<9a T x _ 9x T a _ 
ox ox 
*g^x = ( ^ + A)x ifA ^ T 2Ax 
(written as column vectors) 

(AB) T = B T A T . 

We will now derive each term of the left side of (23). First, let us rewrite L 
in two different ways 

L = ^uj' T Juj' = 2(Gq) T J(Gq) = 2(Gq) T J(Gq) 
and grouping around J 

L = \Q' T JQ' = 2q T (G T JG)q = 2q T (G T JG)q. 



Because J is symmetric, (G T JG) and (G T JG) are also symmetric. So we 
have 

— = 4G T JGq = 2G T J (2Gq) - -2G T Juj' , (24) 

-lD' 

^ = AG T JGq = 2G T J (2Gq) = 2G T Juo' 
oq -—^—y 



and 



±.!*L = 4(2G T J^') - 2G t Juj' + 2G t Jlo'. (25) 
dtdq dV ' K ' 
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1.6.2 Generalized Forces 

A way to find the generalized force F c relative to the coordinates c is to identify 
it in 

SW = F c • Sc. 

(A simple example is the case of a pure translation Sx of a particle, on which 
a force F is applied. The work is then SW = Fg ■ Sx — F ■ Sx. So the generalized 
force Fg is simply F in this case.) 

For a rotation of a rigid body by an angle Sip around an axis n with an 
applied torque T', the work can be written as 

SW=(n-T')5ip \n\ = l- (26) 

This small attitude change can be represented on one side as a small variation 
Sq of the coordinate quaternion q and, on the other side, as a rotation quaternion 
qs operating from the current attitude represented by q (i.e. a composition). 
That is 



q + Sq = q o q s 
\q\ = l \q 6 \ = l 

We do not need to consider the fact that the variation Sq has to preserve the 
norm of q, because it will automatically be satisfied by introducing a constraint 
in the Lagrange formulation. 



On one side we can write 



q + Sq = q o q s 
q o q +q o Sq = q o q oq$ 

(1,0) (1,0) 

qs 

q 5 = (1,0) + qoSq. (27) 

On the other side 

Sip . Sip 
qs = (cos—, sin — n). 

Looking at the imaginary part 



Imj^j = Imjg o Sqj = sin — n ~ — n 
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comparing with (26) 

6W = 

and from (20) 



2 lm{qodq} ■ f ' 



Im{<7 o dq} = GSq 
SW = 2(G8q) ■ f' = 2T' T (G5q) = 2{G T f') T 5q = 2(G T f') -Sq 



F q = 2G T T' 



(28) 



1.6.3 Dynamics 

We have now everything to write the dynamics 



d_dL _dL_ F +x ^£_ 
dt dq dq q dq 



AG T JQ' + 2G T jj' = 2G 1 T' + Aq. 



Trpt 



Left-multiplying by G 



AGG^Juj' + 2GC?Jlo' = 2GC?f' + A Gq 

2W Id Id 

0,'JQ 1 + Jw' = f' 

a' x ju' + jd' = f 
ja' = f ' - a' x jcs. 

This last relation is nothing else than the Euler equation of motion for ro- 
tating body. Together with (18) we obtain the complete dynamics 



q 



J^T'-J-^lj' x JQ 1 ) 
\G T Q'. 



(29) 
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A Derivatives and Quaternions 

A.l Quadratic Form Derivative by a Quaternion 

In order to be able to derive the Lagrangian by the components of q for a 
non-inertial quaternion model, one needs to perform things like 

d^Rw) d(jF R T w) 



<9q <9q 

and also 



d{u T RJR T u) 

But because R — EG T and 



<9q 






( "Si 


9o 


93 


-92 




-ft 


-93 


9o 


9i 




V -93 


92 


-9i 


9o 



E = 



the matrix of the quadratic form to be derived is not constant in q. This 
implies that these operations are no more trivial. However, thanks to the partic- 
ular form of the dependance of R in the components of q, higher order tensors 
can be avoided, as shown in the following. 

A. 1.1 "Single R" Quadratic Form 

By computing the quadratic form and taking the partial derivatives we get 
(placing them in a column vector) 

d^Rw) ( di^FRw) 



<9q V dqi 

/ W! vi q + W! vg qs - w } vs q 2 - w 2 vi qs + w 2 vg qo + w 2 v 3 qi + w s vi q 2 - w 3 v 2 qj + ws vs qo \ 
wi vi qi + wi v 2 q 2 + wi vs qs + w 2 vi q 2 — wg v 2 qi + w 2 vs qo + ws vi qs — ws vg qo — ws vs qi 
-wi V] qg + W! vg qi - wj v 3 qo + w 2 vi qi + wg v 2 q 2 + w 2 v 3 qs + ws vi qo + w s vg qs - w 3 v 3 q 
\ -wi vi qs + W! vg qo + wi v 3 qi - w 2 vi qo - wg vg qs + w 2 v 3 q 2 + w 3 vi qi + w s vg qg + w 3 v 3 qs I 



The vector obtained is quite ugly but one can see that it is linear in q, it 
can thus be rewritten in a matrix- vector product: 



/ vi wi + vg w 2 + vs ws vs w 2 — vg ws 

vs w 2 — vg ws vi wi — vg w 2 — vs ws 
— vs wi + vi ws vi wg + vg wi 

\ vg wi — vi wg vi ws + vs wi 



— vswi + VI ws 

Vt W 2 + Vg W! 
• W 2 - V! W! — Vs Ws 

vg ws + vs w 2 



vg wi — VI wg 

VI W 3 + Vs W! 

vg ws + vs wg 
Vs Ws — vi wi — vg 



w 2 ) 



\ ( QO \ 

11 



A[v, w] 
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By careful inspection of A[v, w], we can identify a structure in the matrix 
that allows a compact notation 



A [v, w] 



w • v (w x v) 

w x v w-tF + vw T — w ■ v I 3 



That is 



djv 1 Rw) 
dq 



= 2A[v,w]c[ 



And because v T R T w = w T Rv we also have 
8{v r R T w) 



dq 



= 2A[w,v\q 



(30) 



(31) 



(32) 



A. 1.2 "Double i?" Quadratic Form 

We are now interested in the derivative of a quadratic form involving RJR T , 
that is, with the q dependent matrix R appearing twice. J is an inertia matrix, 
therefore, J = J T . This time, the vectors on the left ant on the right are the 
same, lets say u. 



ld_ 
2dq 



, dR 
dqi 



u T — JR T u \ . 



Hence 



I A (u T RJR T u) = 2A[u, JR T u]q 

A. 1.3 Properties 

By looking at (30), one may note the following relations 
A[vi + V2,w] = A[v ly w] + A[v 2 ,w] 



(33) 



(34) 



A[v, wi + W2] = A[v, wi] + A[v, w 2 ] 



(35) 



n m 

1=1 3=1 

A[av, (3w\ = a(5A[v, w 



i=l 3 = l 



(36) 
(37) 
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A. 2 Time Derivative of R 

First note that by identification, one can verify that 

G T G = E T E = h- qq T 
with I4 the identity matrix in R 4 . Remember also 

0' = 2GG T = —2GG T with tl'v-- 

and 



Now observe 



lu' = 2Gq = -2Gq. 

n'R T = 2GG T GE T 
= -2GG T GE T 
= -2G(h - qq T )E T 
= -2GE T - 2Gq q T E T 



(Bq) T =0 



We can finally write 



-2GE T = -if. 



if = -n'R T 
it = -m /T = m'. 
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B Speed Composition 



Let be three referentials each designed by 0, 1 and 2. Referential is inertial, 
referential 1 is a rotating one and 2 is a body fixed referential. 
The same vector x can be expressed in any of these referentials; when expressed 
in 0, we will notate it as x°, when expressed in 1 it will be noted x 1 and x 2 in 
referential 2. We will also write x l the quaternion (0, x 1 ). 

Moreover, three quaternions are defined: goi describes relative attitude of refer- 
ential 1 with respect to referential 0, q\ 2 describes relative attitude of referential 
2 with respect to referential 1 and q$ 2 describes relative attitude of referential 
2 with respect to referential 0. 




So we may write 

x° = q (n ox 1 o q al x 1 = q u ° X 2 o q 12 x° = q 02 o x 2 o q 02 
and by substitution 

a;° = qoi ox 1 o q 01 = q Q1 o q 12 o x 2 o q 12 o q 01 = (q 01 o q 12 ) o x 2 o (g 01 o q 12 ) 
we can identify q 02 

902 = 9oi° 9i2- (41) 
Noting ujjj = (0, CSfj) the rotation velocity of the reference frame j relative 
to frame i expressed in the frame j and remembering that Lufj = 2q\j o q^ , we 
may write 
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w 02 = 2?02 ° 902 

= 2(gi2 0501)0 (901 o <7i2 + goi 912) 

= 2gi2 o goi ° 4oi 312 + 2gi2 o q 01 o <7 01 o<j 12 

Id 

= 9i2 o (2g i 9oi ) °5i2 + 2gi 2 o q 12 
. ' v . ' 

<»01 U 12 

= qi2 o Wqi o g 12 + W? 2 

2 2 
= U m + U! 12 . 

That is, we can add consecutive rotation speeds if they are expressed in the 
same referential. 

In the case of the Cubsat, u)q 2 is the satellite's rotation velocity uj' expressed in 
body coordinates in the inertial referential model; we will note it uj' Inertial here. 
On the other hand, uj\ 2 ^ s the satellite's rotation velocity uj' expressed in body 
coordinates in the non-inertial referential model (i.e. in orbital reference frame, 
ORF); we will note it u' NonInertial . 

Wqi is the ORF rotation velocity expressed in the ORF, that is uj , while cjqi 
is the same vector, transformed in the body referential. This transformation is 
performed by R T from the non-inertial model (q~\ 2 in the above developemcnt). 
In other words, we can link the uj' vector from both inertial and non-inertial 
formulations (models) with 

^Inertial = ^Nonlnertial^o + & N on Inertial' (^2) 

This is the speed to be used in computing the kinetic energy for the non- 
inertial model. 
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C Euler Angles to Quaternions 



Three rotations by the Euler angles around each axis can be written as 



cos (tp) — sin (tp) 
sin (-0) cos (tp) 
1 



R$ 



cos (9) sin (0) 

1 
- sin (6>) cos (6) 



Rs = 



1 

cos ((f)) — sin (<fi) 
sin (<j)) cos ((p) 



Combined together, they define the rotation matrix 



R — R<t>ReR%p- 



Those three rotations can also be expressed as quaternion rotations 



308 (§■ 

3in(i. 





cos(i 



sin (| 



q^ 



cos (i tp) 



sin(iV) 



The resulting quaternion can then be obtained by multiplying those three 
together 



20 



cos (i cf) cos (i 9) cos (| V) ~~ si n (| 0) si n (| si n (| V') 
cos (i cos (i 0) sin (i </>) + cos (± 0) sin (| 0) sin (i -;/>) 

cos (i V) cos (\ <A) s i n (| 0) — cos (| ^) sm (\ 4>) si n V') 
cos (i 0) cos (i 0) sin (i -0) + cos (| sin [\ </>) sin (i 0) 

Note: that this result depends on the convention used in the order and choice 
of the Euler angles and rotation axes! [5] 
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